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Direct photons measured by the PHENIX experiment at RHIC 
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Abstract. Results from the PHENIX experiment at RHIC on direct photon production in p + p, d+Au, 
and Au+Au collisions at ^/SNN = 200 GeV are presented. In p + p collisions, direct photon production 
at high pt behaves as expected from perturbative QCD calculations. The p + p measurement serves as a 
baseline for direct photon production in Au+Au collisions. In d+Au collisions, no effects of cold nuclear 
matter are found within the large uncertainty of the measurement. In Au+Au collisions, the production 
of high Pt direct photons scales as expected for particle production in hard scatterings. This supports jet 
quenching models, which attribute the suppression of high pr hadrons to the energy loss of fast partons 
in the medium produced in the collision. Low pr direct photons, measured via e"'"e~ pairs with small 
invariant mass, are possibly related to the production of thermal direct photons. 
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1 Introduction 

Depending on their transverse momentum, direct pho- 
tons convey information about different aspects of ultra- 
relativistic nucleus-nucleus {A + A) collisions. High trans- 
verse momentum (pr) direct photons are produced in early, 
hard parton-parton scatterings by processes like quark- 
gluon Compton scattering (q+g — > q + 7). Unlike scat- 
tered quarks or gluons, direct photons do not interact 
strongly with the medium subsequently produced in the 
collision. Highp^ direct photons (pt^6 GeV/c) therefore 
provide a baseline for measuring medium modifications of 
high Pt hadron production. 

A significant fraction of low px direct photons (1 < 
Pt ^ 3 GeV/ c in central Au+Au collisions at Sjvjv = 200 
GeV) is expected to come from a thermalized medium of 
deconfined quarks and gluons, the Quark-Gluon Plasma 
(QGP), possibly created in A + A collisions [1]. These 
thermal photons carry information about the initial tem- 
perature of the QGP. Thermal photons are produced in 
the QGP as well as in the hadronic gas over the entire 
life time of the collision. The initial temperature can be 
extracted comparing measurements to models which con- 
volute the production rate with the space-time evolution 
of the collision. 

At low and intermediate px {pt < 6 GeV/c in central 
Au+Au collisions at y^s^ = 200 GeV) another significant 
source of direct photons might be the interaction of fast 
partons from jets with thermal partons in the QGP [2] 
like qhard + gqcp — > q + 7. 

Possible modifications of direct photon production from 
cold nuclear matter effects can be measured in c?+Au col- 
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lisions where no medium is created. Direct photon mea- 
surements in p -\- p collisions are a superb test of QCD. 
They constrain the gluon distribution function since the 
gluon is a direct participant of the partonic scattering. 
Furthermore, they provide a baseline for understanding 
direct photon production in A + A collisions. 



2 High Pt direct photons 

2.1 Measurement 

Experimentally, the measurement of direct photons is chal- 
lenging due to a large background from hadron decays like 
7r° — > 77 and ry — > 77. PHENIX measures all those pho- 
ton sources. Photons are measured with the electromag- 
netic calorimeter (EMCal). Neutral pions and 77 mesons 
are measured through their two-photon decay branch [3, 
4] . The EMCal [5] consists of six sectors of a lead scintilla- 
tor and two sectors of a lead glass calorimeter centered at 
midrapidity (77 < 0.35) and covering 1/4 of the azimuthal 
angle. 

The direct photon spectrum is obtained by subtract- 
ing the decay photon spectrum from the inclusive pho- 
ton spectrum. This is gained from photon-like showers 
in the EMCal corrected for contaminations from charged 
hadrons and neutrons. The decay photon spectrum is cal- 
culated from the measured 7r° and 77 spectrum, taking into 
account minor contributions from other hadrons that de- 
cay into photons. 

If the detector occupancy is low as in p + p or d+Au 
collisions, the signal-to-background ratio can be improved 
by event by event tagging of photons that have a matching 
partner as decay photons. 
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Fig. 1. Preliminary direct photon cross section as a function of 
Pt for p+p colhsions at -^s — 200 GeV [6]. The solid curves are 
pQCD predictions [7] for three different scales, which represent 
the uncertainty on the calculation. 

2.2 Results 

2.2.1 p + p collisions 

The preliminary direct photon cross section for p + p col- 
lisions at y/s — 200 GeV is shown in Fig. 1 [6] . It is based 
on an integrated luminosity of 266 nb~^ collected in RHIC 
Run-3. For the whole pr range of 5 to 16 GeV/c the ob- 
served cross section is consistent with a next-to-leading- 
order perturbative-QCD (NLO pQCD) calculation [7]. This 
measurement establishes a reference for direct photon pro- 
duction inA + A. In RHIC Run-5, PHENIX took a p + p 
data set with ten times the Run-3 statistics. 



2.2.2 Au+Au collisions 



PHENIX Au+Au (central collisions): 
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Fig. 2. Nuclear modification factor for direct photon, n , and 
rj production in central Au+Au collisions at y^Srm = 200 GeV. 
For direct photons, a pQCD calculation is used as the reference. 

As shown by PHENIX in RHIC Run-2 [8], at high 
Px direct photon production is consistent with the T^a- 
scaled pQCD expectation also in Au-|-Au collisions. Fig- 
ure 2 compares the nuclear modification factor for direct 
photon, 7r°, and tj production in central Au-|-Au collisions 
at y'Sjvjv = 200 GeV. For direct photons, the pQCD cal- 
culation is used as a reference. The suppression for tt" 
and ?7, supported by the non-suppression of photons, can 
be described with energy loss of partons in the produced 
medium. 

While high pr direct photon production is consistent 
with the pQCD expectation, one should keep in mind that 
there are a number of medium effects that might alter di- 
rect photon production. On the one hand, a significant 
fraction of direct photons is expected to stem from the 
fragmentation of hard-scattered partons into jets. As those 
partons lose energy in the medium, also the fragmenta- 
tion photons should be suppressed. On the other hand, as 
partons lose energy through gluon bremsstrahlung, they 
should also generate photon bremsstrahlung, enhancing 
the direct photon yield. It would be a coincidence if these 
counterbalancing effects exactly canceled. 

Experimentally, contributions from these effects can be 
quantified by studying the azimuthal distribution of direct 
photon production with respect to the reaction plane. If 
fragmentation photons are suppressed, less direct photons 
are expected to be emitted outside of the reaction plane. 
If bremsstrahlung is enhanced, more direct photons will 
be found outside of the reaction plane. 



The high pr direct photon yield in collisions relative 

to p + p is expected to scale with the parton luminosity in 
the overlap region of the two nuclei. The parton luminosity 
is quantified by the nuclear overlap function, Taa- In the 
absence of nuclear effects the nuclear modification factor 



Raa{pt) 



dN/dprl 



A+A 



(Taa)/ X da/dpT\p+p 
is unity for particle production from hard scattering. 



(1) 



2.2.3 d+ Au collisions 

The preliminary direct photon yield for minimum-bias 
d-|-Au collisions at ^Sjvjv = 200 GeV is shown in the up- 
per panel of Fig. 3 along with the yield in p + p collisions. 
The d+Au analysis is based on w 3 billion events sam- 
pled in RHIC Run-3. The lower panel shows the ratio to 
the r^yi-scaled NLO pQCD calculation. The ratio is con- 
sistent with unity over the entire pt range, showing no 
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Fig. 3. Prehminary direct photon invariant yield as a function 
of pt for p + p [6] and minimum-bias d-l- Au coUisions at -^Snn 
= 200 GeV. The solid curves are pQCD predictions [7] for 
three different scales, which represent the uncertainty on the 
calculation. 

indication of cold nuclear matter effects. However, the un- 
certainty of the measurement is large. A high statistics 
d+Au run is planned. 



3 Low Pt direct photons 

3.1 Measurement 

The EMCal measurement at low px suffers from a large 
background from decay photons. Also, the relative energy 
resolution of the EMCal becomes worse at low energy. 
Therefore, a new method [9,10] to measure direct pho- 
tons at low Pt has been employed. It has been carried out 
in heavy ion experiments for the first time. The method 
is based on the measurement of pairs of e"'"e^, which are 
identified with the PHENIX Ring Imaging Cherenkov De- 
tector (RICH). The basic idea is that any source of real 
photons also produces virtual photons that decay into 
e~^e~ with small invariant mass. An example is the tt'^ 
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Fig. 4. Invariant-mass distribution of virtual photons from the 
tt" and i] Dalitz decay as well as from direct photons according 
to Eq. 3. It is illustrated how the various contributions decrease 
to a fraction R when going to higher invariant mass with the 
tt" contribution exhausting. 

Dalitz decay (tt^ — > 7e+e~). This internal conversion 
method is based on two assumptions: first, the ratio of 
direct-to-all photons is the same for real and virtual pho- 
tons at small invariant mass close to zero (m-y < 30MeV) : 
^direct/^inci. = Idirect / Unci. ■ Sccoud, the mass distribu- 
tion can be described by the KroU-Wada formula [11], 
which has been established to describe the Dalitz decay ^ : 



2^(1 + I Pi'^ee) I (1 - irf ) , 

The mass distribution according to Eq. 3 is depicted 
in Fig. 4. For e^e^ pairs from tt*^ and rj Dalitz decays the 
yield is suppressed towards higher niee due to the mass, 
M, of the parent meson in the phase space factor (1 — 
ml^/M'^Y whereas no suppression takes place for e+e" 
pairs from virtual direct photons as long as mee vW ■ 
For the small invariant masses considered here the form 
factor is assumed to be unity in all cases. 

The key advantage of this method is the greatly im- 
proved signal-to-background ratio, which is achieved by 
eliminating the contribution from tt*^ Dalitz decays when 
the invariant mass is increased. In addition, the electron 
measurement through charged particle tracking provides 
a better energy resolution at low pT than the EMCal pho- 
ton measurement. This excellent energy resolution at low 
Pt combined with little material upstream of the detector, 
where photons could convert generating background e^e~ 
pairs, makes this measurement feasible in PHENIX. 

^ Since the virtual photons decay in the medium, this re- 
lation might be slightly modified. This would not affect the 
significance of the observed excess of direct photons, only its 
translation into an absolute yield of real direct photons. 
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The experimentally observed quantity is the ratio of 
e~^e~ pairs in an invariant mass bin where the 7r° Dalitz 
decay is largely suppressed to the yield at an invariant 
mass close to zero: i?data = Nf ^-^^"^"^ /N<^°^''^ . If there 
is no direct photon signal, i?data = -Rhldron' -'^data 
can be calculated^ from Eq. 3 based on the known ratio 
rj/n^ = 0.45 ± 0.1. An excess i?data > -Rhadron translates 
into the fraction of virtual direct photons at close to zero 
invariant mass according to: 



7d 



irect 



p pcalc 

-"■data -"-hadron 



^incl. "'Mirect7 '"'hadron 



pcalc pcalc 

direct7 -"-t..-!- 

This can be derived starting from 

D Ar90-300MeV /ArOOMeV 



'/K 



(3) 



(4) 



with 



-^ee — -^hadron H" -^direct 7 



(5) 

Dcalc jyOOMeV , ocalc lyOOMeV 

p hadron hadron direct 7 direct 7 

-n-data — jyOOMeV , jyOOMeV ^ ) 

hadron ' direct 7 

jyOOMeV 

_ pcalc I /-pcalc _ pcalc n direct 7 

-""hadron "r V-'^'direct 7 -"'hadron/ ^<;30MeV ^ ' 

The direct photon spectrum is then obtained by multiply- 
ing l^irect/^inci ^y the iuclusive photon spectrum mea- 
sured with the EMCal. 

The full 2004 data set of about 900 M minimum bias 
events is analyzed. Events and centrality are selected as 
described in [12]. Electrons in the central arms are iden- 
tified by matching charged particle tracks to clusters in 
the EMCal and to rings in the RICH detector. To ob- 
tain a clean invariant-mass distribution of e^e~ pairs, 
pairs originating from photon conversions in the beam 
pipe or detector material are rejected based on their orien- 
tation with respect to the magnetic field. The combinato- 
rial background is removed by an event-mixing technique. 
The uncertainty of the ry-to-Tr^ ratio of about 20 % [13] is 
the main source of uncertainty, translating into an uncer- 
tainty of 20 % of the measured direct photon yield. Other 
sources are the EMCal-measured inclusive photon yield 
(10 %) and the e+e~-pair acceptance (5 %). The total 
systematic uncertainty is 25 %. 



based on a subset of the RHIC Run-4 data set [14]. ^ 
While the internal conversion measurement results in a 
significant direct photon signal of about 10 % above the 
decay photon backgromid for 1 < < 5 GeV/c, the 
EMCal measurement does not yield a significant direct 
photon signal below = 3 GeV/c, but agrees with the 
internal conversion measurement within the uncertainties 
over the entire range. 

Figure 5 b) shows the direct photon invariant yield 
from the internal conversion measurement [14] and com- 
pares it to theoretical calculations. With large significance, 
a direct photon spectrum is obtained for I < px < 5 
GeV/c. The spectrum lies significantly above a TA^-scaled 
NLO pQCD calculation [7] for pr <3 GeV/c. The pQCD 
calculation indicates the contribution from hard scatter- 
ings. The lines show the scale uncertainty of the calcula- 
tion. However, it is not clear how meaningful the compar- 
ison is down to this low pt where the pQCD calculation 
reaches its limit of applicability. It is planned to replace 
the pQCD calculation by a reference measurement of di- 
rect photon production in p + p collisions with the same 
method as in Au-|-Au. 

The excess of the measured direct photon spectrum 
above the pQCD calculation can be described by models 
that allow for a significant contribution of thermal pho- 
tons. In order to describe thermal photon production in 
A + A collisions, the entire space-time evolution has to be 
accounted for including the hadronic phase. This is done 
in hydrodynamical models, which assume local thermal 
equilibrium. An important free parameter in such mod- 
els is the initial temperature of the fireball. Figure 5 b) 
compares the measured direct photon spectrum to a 2-|-l 
hydrodynamical model [15] for thermal-photon emission 
with an average initial temperature of Tq^'^ = 360 MeV 
(•ymax ^ MeV) and a formation time of tq = 0.15 
fm/c. The model underpredicts the data for pT >3 GeV/c. 
The data can be described when both sources, thermal and 
pQCD, are combined. Various other calculations come to 
similar descriptions with initial temperatures in the range 
370 < Ti < 570 MeV [2,16,17]. The temperatures are 
significantly above the critical temperature for the QGP 
phase transition of Tc = 170 MeV. However, the obtained 
temperature is only meaningful if the excess above pQCD 
is of thermal origin. Besides hard photons, also jet-plasma 
interactions may play a significant role [2]. 



3.2 Results 

Figure 5 a) shows the direct photon signal for central 
Au-|-Au collisions at ^/s]^ = 200 GeV in terms of the 
double ratio (7/7r°|meas.)/(7/7r°|bckgrd.)- This ratio indi- 
cates a direct photon excess above the decay photon con- 
tribution as an enhancement above 1. In this figure, two 
measurements are compared: the result from the internal 
conversion method as described above [14] and a prelim- 
inary analysis of the conventional EMCal measurement 

^ Besides tt" and r), decay photons from a cocktail of hadrons 
are considered here. 



4 Conclusions 

The direct photon cross section in p + p collisions at ^/s 
= 200 GeV is found to be consistent with a NLO pQCD 
calculation over the entire range of the measurement from 

5 GeV < Pt < 16 GeV/c. In d+Au collisions, no indica- 
tions for cold nuclear matter effects are found within the 

^ The preliminary result from the measurement of vir- 
tual photons is presented in the same figure in terms of 
^direct /"find. + 1, which also iudicates a direct photon excess 
as an enhancement above 1. Note that the two quantities are 
not exactly equivalent. 
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Fig. 5. a) Direct photon excess for the conventional and the internal-conversion measurement in central Au+Au collisions [14]. 
b) Direct photon spectrum from the latter [14] compared to pQCD [7], thermal-photon [15], and the sum of both calculations. 



large uncertainties of the measurement. A high statistics 
d-l-Au run is planned. Highp^ direct photon production in 
Au-|-Au collisions is found to scale with the nuclear overlap 
function Taa, establishing medium effects as the cause for 
the hadron suppression observed in central Au-|-Au colli- 
sions. 

At low pt (1 < Pt < 5 GeV/c) a significant direct 
photon excess above the decay photon background is ob- 
served employing a new method in heavy ion collisions, 
which measures low invariant mass e'^e~ pairs from di- 
rect photon internal conversions. The signal appears to 
be above pQCD calculations for < 3 GeV/c, whose 
applicability is questionable, however, at this low px- If a 
thermal photon source modeled by hydrodynamical calcu- 
lations is added to the pQCD calculation, the data can be 
described over the entire pt range. Other direct photon 
sources like jet-plasma interactions might play a signifi- 
cant role as well. 
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